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Abstract 

We study geometric consistency relations between angles on S-dimensional (3D) circu- 
lar quadrilateral lattices — lattices whose faces are planar quadrilaterals inscribable into 
a circle. We show that these relations generate canonical transformations of a remarkable 
"ultra-local" Poisson bracket algebra defined on discrete 2D surfaces consisting of circu- 
lar quadrilaterals. Quantization of this structure leads to new solutions of the tetrahedron 
equation (the 3D analog of the Yang-Baxter equation). These solutions generate an infinite 
number of non-trivial solutions of the Yang-Baxter equation and also define integrable 3D 
models of statistical mechanics and quantum field theory. The latter can be thought of 
as describing quantum fluctuations of lattice geometry. The classical geometry of the 3D 
circular lattices arises as a stationary configuration giving the leading contribution to the 
partition function in the quasi-classical limit. 
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1 Introduction 



Currently the quantum integrability is mostly understood as a purely algebraic phenomenon. 
It stems from the Yang-Baxter equation [1,2] and other algebraic structures such as the affine 
quantum groups [3,4] (also called the quantized Kac-Moody algebras), the Virasoro algebra [5] 
and their representation theory. It is, therefore, quite interesting to learn that some integrable 
models of statistical mechanics and quantum field theory arise also in quantization of models of 
discrete differential geometry. Recently, we have shown [6, 7] that the two-dimensional integrable 
model associated with the Faddeev-Volkov solution [8-10] of the Yang-Baxter equation describes 
quantum fluctuations of circle patterns [11] connected with the Thurston's discrete analogues of 
the conformal transformations of the two-dimensional plane [12]. The Faddeev-Volkov model is 
an Ising-type model with continuous spin variables, which are interpreted as fluctuating radii of 
the circles. It contains a free parameter, which can be identified with the Planck constant (in 
the language of Euclidean quantum field theory) or with the temperature (in the language of 
statistical mechanics). The classical geometry of the circle patterns is described by stationary 
configurations giving the leading contribution to the partition function in the quasi-classical (or 
the zero-temperature) limit. 

In this paper we consider similar connections between integrable three-dimensional (3D) 
quantum systems and integrable classical models of 3D discrete differential geometry. The analog 
of the Yang-Baxter equation for integrable quantum systems in 3D is called the tetrahedron 
equation. It was introduced by Zamolodchikov in [13, 14] (see also [15-21] for further results in 
this field, used in this paper). Similarly to the Yang-Baxter equation the tetrahedron equation 
provides local integrability conditions which are not related to the size of the lattice. Therefore 
the same solution of the tetrahedron equation defines different integrable models on lattices 
of different size, e.g., for finite periodic cubic lattices. Obviously, any such three-dimensional 
model can be viewed as a two-dimensional integrable model on a square lattice, where the 
additional third dimension is treated as an internal degree of freedom. Therefore every solution 
of the tetrahedron equation provides an infinite sequence of integrable 2D models differing by 
the size of this "hidden third dimension". Then a natural question arises whether known 2D 
integrable models can be obtained in this way. A complete answer to this question is yet 
unknown. So far only two different (but related) examples of such correspondence have been 
constructed. The first example, connected with the generalized chiral Potts model [22,23], was 
found in [16]. The second example was recently found in [21]. The corresponding solution of the 
tetrahedron equation constructed in [21] reproduces all two-dimensional solvable models related 
to finite-dimensional highest weight representations for all quantized affine algebras Uq{sln), 
n = 2,3, ... , oo, where the rank n coincides with the size of the hidden dimension. Here we 
unravel yet another remarkable property of the same solution of the tetrahedron equation. We 
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show that it can be obtained from quantization of geometric integrabihty conditions for the 3D 
circular lattices — lattices whose faces are planar quadrilaterals inscribable into a circle. 

The 3D circular lattices were introduced [24] as a discretization of orthogonal coordinate 
systems, originating from classical works of Lame [25] and Darboux [26]. In the continuous 
case such coordinate systems are described by integrable partial differential equations (they 
are connected with the classical soliton theory [27,28]). Likewise, the quadrilateral and circular 
lattices are described by integrable difference equations. The key idea of the geometric approach 
[24,29-36] to integrability of discrete classical systems is to utilize various consistency conditions 
[37] arising from geometric relations between elements of the lattice. It is quite remarkable that 
these conditions ultimately reduce to certain incidence theorems of elementary geometry. For 
instance, the integrability conditions for the quadrilateral lattices merely reflect the fact of 
existence of the 4D Euclidean cube [30]. Here we present these conditions algebraically in a 
standard form of the functional tetrahedron equation [20]. The latter serves as the classical 
analog of the quantum tetrahedron equation, discussed above, and provides a connecting link 
to integrable quantum systems. 

We study relations between edge angles on the 3D circular quadrilateral lattices and show 
that these relations describe symplectic transformations of a remarkable "ultra-local" Poisson 
algebra on quadrilateral surfaces (see Eq. (f25l) ). In Section [3] we formulate a variational principle 
which generates these angle relations and explicitly calculate the Lagrangian form of the action 
functional 5*^'^'^ (curiously, it expressed through the Lobachevsky function). Next, in Section U] 
we quantize this structure and obtain two solutions of the tetrahedron equation. One of them 
was previously known [21], but another one, given in Section [4.31 and [4.41 is new. The solutions 
are used to define two different integrable models of statistical mechanics and quantum field 
theory on the cubic lattice. Their partition functions depend on the quantum parameter h (the 
Planck constant). For both of these models the quasi-classical limit of the partition function 



is governed by the above classical action S^'^^\ determined by the angle geometry of the circular 
lattices (more precisely, the two models lead to two different regimes of the classical action 
connected by an analytic continuation). In this paper we only state our main results, the details 
will be presented elsewhere [38]. 

2 Discrete differential geometry: "Existence as integrability" 

In this section we consider classical discrete integrable systems associated with the quadrilateral 
lattices. There are several ways to extract algebraic integrable systems from the geometry 
of these lattices. One approach, developed in [30,32,39-41], leads to discrete analogs of the 
Kadomtsev-Petviashvili integrable hierarchy. Here we present a different approach exploiting 
the angle geometry of the 3D quadrilateral lattices. 

2.1 Quadrilateral lattices 

Consider three-dimensional lattices, obtained by embeddings of the integer cubic lattice Z'^ into 
the A^-dimensional Euclidean space M^, with A*" > 3. Let x{m) G M^, denote coordinates of the 
lattice vertices, labeled by the 3-dimensional integer vector m = miei+m2e2+Tn^es G Z^, where 
ei = (1,0, 0), 62 = (0,1,0) and 63 = (0,0,1). Further, for any given lattice vertex xq = x{m), 
the symbols Xi = x{m + Cj), Xij = x{m + ei + ej), etc., will denote neighboring lattice vertices. 




h 







(1) 
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The lattice is called quadrilateral if all its faces {xQ,Xi,Xj,Xij) are planar quadrilaterals. The 
existence of these lattices is based on the following elementary geometry fact (see Fig. 1) [30], 



3=12 




Figure 1: An elementary hexahedron of a cubic quadrilateral lattice. 

Consider four points Xq,Xi,X2,X3 in general position in R^, > 3. On each of the three 
planes {xo,Xi,Xj), 1 < i < j < 3 choose an extra point Xij not lying on the lines {xQ,Xi), 
{xQ,Xj) and {xi,Xj). Then there exist a unique point X123 which simultaneously belongs to the 
three planes (a^i, 0^12, a^is), (2:2, a;i2, a;23) and {x3,xi3,X23)- 

The six planes, referred to above, obviously lie in the same 3D subspace of the target space. 
They define a hexahedron with quadrilateral faces, shown in Fig. [TJ It has the topology of the 
cube, so we will call it "cube" , for brevity. Let us study elementary geometry relations among 
the angles of this cube. Denote the angles between the edges as in Fig. [2l Altogether we have 




Figure 2: The "front" (a) and "back" (b) faces of the cube in Fig. [T] and their angles. 



6 X 4 = 24 angles, connected by six linear relations 

aj + Pj + 7j + 6j = 27r, a'j + + + dj = 2tt, j = 1, 2, 3, (2) 

which can be immediately solved for all "5's". This leaves 18 angles, but only nine of them are 
independent. Indeed, a mutual arrangement (up to an overall rotation) of unit normal vectors 
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to six planes in the 3D-space is determined by nine angles only. Once this arrangement is fixed 
all other angles can be calculated. Thus the nine independent angles of the three "front" faces 
of the cube, shown in Figl2^, completely determine the angles on the three "back" faces, shown 
in FigllJ), and vise versa. So the geometry of our cube provides an invertible map for three 
triples of independent variables 

^123: {aj,Pj,lj} ^ {aj,f3j,lj}, j = 1,2,3. (3) 

Suppose now that all angles are known. To completely define the cube one also needs to specify 
lengths of its three edges. All the remaining edges can be then determined from simple linear 
relations. Indeed, the four sides of every quadrilateral are constrained by two relations, which 
can be conveniently presented in the matrix form 

V C{A) D{A) J \ su^ I 

\ sm<3 sino / 

where A = {a,/3, 7,5} denotes the set of angles and ip,£g,ip/,iqi denote the edge lengths, 
arranged as in Figl3j Note that due to ([2]) the entries of the two by two matrix in (jj]) satisfy 
the relation 

AD - BC = {AB -CD)/{DB - AC). (5) 

Assume that the lengths (p,iq, £r, on one side of the two pictures in FigH] are given. Let 
us find the other three lengths ip',iq', Ir' on their opposite side, by iterating the relation Q. 
Obviously, this can be done in two different ways: either using the front three faces, or the back 
ones — the results must be the same. This is exactly where the geometry gets into play. The 
results must be consistent due to the very existence of the cube in Fig. [1] as a geometric body. 
However, they will be consistent only if all geometric relations between the two sets of angles 






X23 



Figure 3: The angles A — {a, (3, 7, S} and sides ip, £q, £pi , £q> of a quadrilateral and the oriented rapidity 

lines. 



in the front and back faces of the cube are taken into account. To write these relations in a 
convenient form we need to introduce additional notations. Note, that Figl3] shows two thin 
lines, labeled by the symbols "p" and "g" . Each line crosses a pairs of opposite edges, which we 
call "corresponding" (in the sense that they correspond to the same thin line). Eq.([3]) relates 
the lengths {ip, Iq) of two adjacent edges with the corresponding lengths {I'p, I'q) on the opposite 
side of the quadrilateral. 

Consider now FigH^ which contains three directed thin lines connecting corresponding edges 
of the three quadrilateral faces. By the analogy with the 2D Yang-Baxter equation, where similar 
arrangements occur, we call them "rapidity" line^. We will now apply ([!]) three times starting 

^However, at the moment we do not assume any further meaning for these hnes apart from using them as a 
convenient way of labeling to the corresponding (opposite) edges of quadrilaterals. 



5 




from the top face and moving against the directions of the arrows. Introduce the following three 
by three matrices 

/ A B 0\ / A B \ /100\ 

XpgiA) = i C D \ , XpriA) = 1 , X,r{A) = I A B \ , (6) 
\ I J \ C D J \0 C D J 

where A,B,C,D are defined in ([ID and their dependence on the angles A = {a, (3,^,5} is 
implicitly understood. It follows that 

{i'pj'.j'rY = Xp,{Al) Xpr{A2) XgriAs) {ipjgjrf (7) 

where 

Aj = {a,,P,,j,,5j}, J = 1,2,3, (8) 

the lengths ip,iq, . . . are defined as in FigHl and the superscript "t" denotes the matrix trans- 
position. Performing similar calculations for the back faces in FigUja and equating the resulting 
three by three matrices, one obtains 

Xpq{Al) Xpr{A2) XqriA-s) = Xqr{A'^) Xpr{A2) Xpq{A'i) . (9) 

where 

A'j = {a'j,P'j,-f'j,6'^}, j = l,2,3. (10) 

This matrix relation contains exactly nine scalar equations where the LHS only depends on the 
front angles dH), while the RHS only depends on the back angles (fTOj) . Solving these equations 
one can obtain explicit form of the map ([3]). The resulting expressions are rather complicated 
and not particularly useful. However the mere fact that the map ([3]) satisfy a very special Eq. ([9|) 
is extremely important. Indeed, rewrite this equation as 

XpqiAl) Xpr{A2) XqriAs) = 7^l23 [XqriAs) Xp,(^2) Xpq{Al)) (11) 

where TZi23 is an operator acting as the substitution ([3]) for any function F(^i,^2,^3) of the 
angles, 

ni23[F{AuA2,A3)) = F{A[,A'2,A's) (12) 
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Then, following the arguments of [19], one can show that the map ^ satisfies the functional 
tetrahedron equation [20] 

"7^123 ■ '^145 • ^246 ' ^356 = ^356 ' ^246 ' ^145 ' '^123 ) (13) 

where both sides are compositions of the maps ([3|), involving six different sets of angles. Alge- 
braically, this equation arises as an associativity condition for the cubic algebra To discuss 
its geometric meaning we need to introduce discrete evolution systems associated with the map 

2.2 Discrete evolution systems: "Existence as integrability" 

Consider a sub-lattice L of the 3D quadrilateral lattice, which only includes points x{m) with 
mi,m2,rn-3 > 0. The boundary of this sub-lattice is a 2D discrete surface formed by quadrilat- 
erals with the vertices x{m) having at least one of their integer coordinates mi, ?Ti2, ma equal to 
zero and the other two non-negative. Assume that all quadrilateral angles on this surface are 
known, and consider them as initial data. Then repeatedly applying the map ([3|) one can calcu- 
late angles on all faces of the sub-lattice L, defined above (one has to start from the corner a;(0)). 
The process can be visualized as an evolution of the initial data surface where every transforma- 
tion ([3]) corresponds to a "flip" between the front and back faces (Fig. [2]) of some cube adjacent 
to the surface. This makes the surface looking as a 3D "staircase" (or a pile of cubes) in the 
intersection corner of the three coordinate planes, see Fig. [5] showing two stages of this process. 
Note, that the corresponding evolution equations can be written in a covariant form for an arbi- 




Figure 5: Visualization of the 3D "staircase" evolution. 

trary lattice cube (see Eq. ()26p below for an example). It is also useful to have in mind that the 
above evolution can be defined purely geometrically as a ruler- and- compass type construction. 
Indeed the construction of the point X123 in Fig. [1] from the points a^o, a^i, 3^2, 333, a3i2, 3^13, iC23 
(and that is what is necessary for flipping a cube) only requires a 2D-ruler which allows to draw 
planes through any three non-collinear points in the Euclidean space. 

Similar evolution systems can be defined for other quadrilateral lattices instead of the 3D 
cubic lattice considered above. Since the evolution is local (only one cube is flipped at a time) 
one could consider finite lattices as well. For example, consider six adjacent quadrilateral faces 
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covering the front surface of the rhombic dodecahedrorH shown in Fig. [UJ Suppose that ah angles 
on these faces are given and consider them as initial data. Now apply a sequence of four maps 
([3]) and calculate angles on the back surface of the rhombic dodecahedron. This can be done in 
two alternative ways, corresponding to the two different dissection of the rhombic dodecahedron 
into four cubes shown in FigEl The functional tetrahedron equation (fT3l) states that the results 



Figure 6: Two dissections of the rhombic dodecahedron into four quadrilateral hexahedra. 

will be the same. Thereby it gives an algebraic proof for the equivalence of two "ruler-and- 
compass" type constructions of the back surface of the dodecahedron in Fig. [6l Can we also 
prove this equivalence geometrically? Although from the first sight this does not look trivial, 
it could be easily done from the point of view of the 4D geometry. The required statement 
follows just from the fact of existence of the quadrilateral lattice with the topology of the 4D 
cube [30]. The latter is defined by eight intersecting 3-planes in a general position in the 4-space. 
The two rhombic dodecahedra shown in Fig. [6] are obtained by a dissection of the 3-surface of 
this 4-cube, along its 2-faces, so these dodecahedra must have exactly the same quadrilateral 
2-surface. Thus the functional tetrahedron equation (jl3p . which plays the role of integrability 
condition for the discrete evolution system associated with the map ([3]), simply follows from the 
mere fact of existence of the 4-cube, which is the simplest 4D quadrilateral lattice. For a further 
discussion of a relationship between the geometric consistency and integrability see [37]. 

Note that, to our knowledge, the linear problem ^ for the lengths of the quadrilateral with 
coefficients depending on the angles was not hitherto considered. Much attention was given to 
the linear problem for coordinates of lattice vertices introduced in [30]. In the notations of Fig. [3] 
it reads 

X23-X0 = Y23{X2-X0) + Y,2{xs-Xo), ¥23 = '^^, Ys2 = '-^^j , (14) 

sm a ip sm a iq 

where Xq,X2,X3,X23 are 3D coordinate of the vertices. The evolution equations [30,40] for the 
coefficients Yij mix lengths and angles in a complicated way . To the contrary, our map ([3D only 
involve angle variables. It describes internal geometry of the quadrilateral lattice independently 
of its embedding into the target space (the latter, of course, depends on the length variables). 

Finally, note that the same map ([3]) also arises if one considers a linear problem for unit 
vectors associated with directions of the lattice edges. These vectors satisfy the relation (in the 
notation of Figl3|) 

<-'>'(::) <^^' 



^It is worth noting that the most general rhombic dodecahedron with quadrilateral faces can only be embedded 
into (at least) the 4D Euclidean space. 
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where X is the two by two matrix defined in ^ and 



_ X2, - Xq _ Xq - X2 I _ X23 - X2 , _ X^ - X23 , , 

— p ^ — p 5 Vp — , — . (ibj 

lip f.q r.p f-q 

The property dS} follows from the condition = 1. 



3 Variational principle for the 3D circular lattices 
3.1 Poisson structure of circular lattices 

The 3D circular lattice [24,33,34] is a special 3D quadrilateral lattice where all faces are circular 
quadrilaterals (i.e., quadrilaterals which can be inscribed into a circle). The existence of these 
lattices is established by the following beautiful geometry theorem due to Miquel [42] (see Fig. [Tj) 




Figure 7: Miquel configuration of circles in 3D space, an elementary hexahedron and its circumsphere. 

Miquel theorem. Consider four points Xo,Xi,X2,x-3 in general position in M^, > 3. On 
each of the three circles c{xQ,Xi,Xj), 1 < i < j < 3 choose an additional new point Xij. Then 
there exist a unique point X123 which simultaneously belongs to the three circles c{xi,xi2,xi3), 

c{x2,Xi2,X2-i) and C{X3,X13,X23)- 

It is easy to see that the above six circles lie on the same sphere. It follows then that every 
elementary "cube" on a circular lattice (whose vertices are at the circle intersection points) is 
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inscribable into a sphere, see Fig. [71 The general formulae of the previous subsection can be 
readily specialized for the circular lattices. A circular quadrilateral has only two independent 
angles. In the notation of Fig. [3] one has 



7 = vr - /3, 



TT — a . 



(17) 



Due to the Miquel theorem we can simply impose these restrictions on all faces of the lattice 
without running to any contradictions. The two by two matrix in takes the form 



X 




det X = 1 



(18) 



where we have introduced new variables 

k = (csca)sin/3, a = (esc a) sin(a + /?) , a* = (esc a) sin(a — /3), (19) 
instead of the two angles {a, (3}. Note that the new variables are constrained by the relation 

aa* = l-k'^ . (20) 

Conversely, one has 



a — a 



cos a 



cos P 



a + a* 



(21) 



2A; ' 2 

Let the variables {kj,aj,aj}, {kj,a'j,aj'}, j = 1,2,3, correspond to the front and back faces of 
the cube. The map ([3]) then read explicitly 



n 



123 



k^a^ ek^a2a^, 



(k^alY 
(02)' = alal + ek^k^a2 



(02)' = a^ag + £^1/0302, 



(22) 



■2' 



where e = +1 and 



k2 = yr^a^ 



(23) 



At this point we note that exactly the same map together with the corresponding equations 
([9]) and (jl3p were previously obtained in [21]. Moreover, it was discovered that this map is a 
canonical transformation preserving the Poisson algebra 



{fflj,ci. } — 2 6ij k^ , {ki,aj} — SijkiCii , {A;j,a } 



6ij ki Clj^ , 



i,j = 1,2,3 



(24) 



where kf = 1 — aia*. Note that variables k,a,a* on different quadrilaterals are in involution. 
The same Poisson algebra in terms of angle variables reads 



{ai,l3j} = 5ij , {ai,aj} = {(3i,(3j} = 



(25) 



This "ultra-local" symplectic structure trivially extends to any circular quad-surface of initial 
data, discussed above. To resolve an apparent ambiguity in naming of the angles, this surface 
must be equipped with oriented rapidity lines, similar to those in Fig. [^. In addition, the angles 



■^We refer the reader to our previous paper [6] where the relationship between the rapidity graphs and quadri- 
lateral lattices is thoroughly discussed. 
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for each quadrilateral should be arranged as in Fig. [31 Then one can assume that the indices 
i,j in (j25]) refer to all quadrilaterals on this surface. 

Thus, the evolution defined by the map ()22p is a symplectic transformation. The correspond- 
ing equations of motion for the whole lattice (the analog of the Hamilton-Jacobi equations) can 
be written in a "covariant" form. For every cube define 

^32 = ai, A23 = al, ^31 = 4, Ai3 = a2, ^21 = 03, ^12 = ^3 (26) 

where Ajk stands for Ajk{m), where m is such that x(m) coincides with the coordinates of the 
top front corner of the cube (vertex xq in Fig{T]). Let be the shift operator Tfcj4jj(m) = 
Aij{m + ek)- Then 

TkAij = — - , Kij = Kji = y^l — AijAji , (27) 

where (i, j, k) is an arbitrary permutation of (1, 2, 3) and 

fi=ri, f2 = T^\ T^ = T^- (28) 

Note that Eq.(l271) also imply 

{TkKij)Kkj = {TiKkj)Kij . (29) 

Remarks. The equations (j27p have been previously obtained in [34], see Eq.(7.20) therein. 
The quantites Aij in (j27p should be identified with the rotation coefficients denoted as (3ij in [34] . 
The same equations ([27]) are discussed in §2.2 of [37], where one can also find a detailed bibli- 
ography on the circular lattices (we are indebted to A.I.Bobenko for these important remarks). 

3.2 Variational principle: the Lagrangian and action. 

Every canonical transformation 

Pi, Qi ^ Pi, Q\ : ^dPiN dQi = Y^dPlA dQ'i (30) 

i i 

defines a generating function C{Q, Q'): 

^' = aiS ■ ^ = -«? ' <!£ = E (^'.''ft - i^'^QS) . (31) 

For discrete (in time) canonical transformations, this generating function coincides with the 
Lagrangian density of the action functional. There are many ways to choose canonical variables; 
the resulting Lagrangians differ by equivalence transformations. A convenient choice is the 
variables log ki and ^ log Vi , 

Vi = ^, { logki , logvj} =26ij . (32) 

The Lagrangian density C = C{y,v'), associated with one cube, is defined as 

3 

2 



1 ^ 

dC{v,v') = - ^ ( log kj d log - log k'^ d log v'^ ) , (33) 
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where v = (fi,f2,f3) and similarly for v'. The equations of motion (|22|) imply the relations 



V ''l"3 



1-^ 



(34) 



Introduce new variables 



1 1 1 V 

"i^a/ \ "i-'a/V "i"3 







(- 


4 


"1^3/ 


v^y 




^) 


(- 




"i"3y 


v^y 



,/ 



-2^2 ^ ""2 -2^1 ^ 2no ^ ^2 -2^3 ^ ^2 



(35) 



2Q'2 ^ ^2 2n[ ^ ^2 2QJ, ^ "^2 2% ^ ''^2 

fit's f^fS Vj^Wg flVg 



O'fc = Ofc - (17o + ^^i +^^2 + f^3)/2, fc = 0,l,2,3. (36) 



such that 



Remarkably, the differential (j33p depends on four independent variables only. Indeed rewriting 
dC{v,v') in terms of six independent variables vi,v^, one can easily see that all 

its dependence on vi, V3 drops out, 

1 / sinh^ sinh^ ^. \ 

= - log sinhO- sinh^sinhO- ) '^^^ ' ^'^^ 

Integrating the last equation, one obtains 

3 

2 



Civ,v') = l^[A{nk)+A{n',)) , (38) 

fc=0 

where A(x) is the modified Lobachevsky function, 

A{x)= I log|2sinhy|dy . (39) 
Jo 

The variables $1^,$!^ in the RHS of (j38p are understood as functions of v,v', defined in (j35p . 

In the above derivations of the equations of motion (j34p and the Lagrangian density ()38p . 
we considered an isolated cube and the associated face variables {kj,Vj}, {kj,Vj}, j = 1,2,3. 
Restoring now the coordinates m € for the whole lattice (see the definitions before Fig{T]), 

Vj vAm), f' vAm + Cj), 

' , , \ j = 1,2,3 (40) 

kj kj(m), kj kj[m + Cj) 

one has for the total action 

5(^0({^}) = ^ . (41) 
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where Vj{m) = Vj{m + ej). The variational principle for this action leads to the Lagrangian 
equations of motion with respect to the variables Vj{m). To write them in a covariant form it is 
convenient to define new variables 

ui{m) = vi{m) , U2{rn) = l/v2{m) , U3{m) = 1^3 (m) . (42) 

For an arbitrary lattice site m G denote 

uf^ = Ui{m), u'p =Ui{m + ej), uf^'^ =Ui{m + ej + ek), (43) 

where the indices i,j, k independently take any of the three values 1, 2, 3. Then the Lagrangian 
equations of motion, determining the stationary point of the action ()4ip take the form 

(1 - u\ 'u)"ul') (1 - u\ 'u) 'ul ') K 'u^'ul ') (1 - u\ 'u) 'ul ') ^^^^ 



where (i, j, k) is any permutation of (1, 2, 3). Note that similar, but different, equations of motion 
arose previously in [34,43]. 

It is sometimes convenient to use the /c- variables, kj{m), instead of the v-variables. The 
corresponding Lagrangian density for one cube is obtained from (|38|) with a Legendre transform 

_ 1 3 3 

dk, k') = --Y, log log V, + C{v, v') + -Yl log k'^ log . (45) 
i=i j=i 

where the variables v^v' are now considered as functions of k^k' implicitly defined by (j34p . The 
total action 

= ^ £(A;(m),A;'(m)) . (46) 

obviously coincides with (|4ip . up to boundary terms. Note that the Lagrangian C{v,v') has a 
gauge symmetry: it depends only on four combinations of the six independent variables Vj,Vj 
(see (|35|36p above). As a result the variables kj, k'- in C{k, k') are constrained by two relations 

which correspond to (j29p . Here we prefer to avoid Lagrangian multipliers and deal with the 
restricted configuration space directly. The local constraint (j47p . taken for all lattice sites, leads 
to certain restrictions for boundary values of /c's and makes the Lagrangian equations of motion 
more complicated than those in the i^-variables. For instance, the local equations of motion 
in the fc- variables involve eight adjacent hexahedra, while those in w-variables, given by (I44p . 
involve only two. Nevertheless, the Dirichlet problem for the action (j46p is well posed [38]. 
Note also that the logarithmic terms in ()45p can be written as 

33 3 

- log kj log Vj + log k'j log v'j = -Y {^3 log(2 sinh Q^) + Q'^ log(2 sinh Q'^)^ . (48) 
j=i j=i j=o 
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4 Quantization. 



4.1 Tetrahedron equation 

In this section we construct two different, but related quantizations of the map (j22p and obtain 
two solutions of the full quantum tetrahedron equation (see Eq. ()56p below). In both cases we 
start with the canonical quantization of the Poisson algebra ()25p . 



[ai,l3j]=^h6ij, [a^,aj]=0, /J^-] = , (49) 

where h is the quantum parameter (the Planck constant) and ^ is a numerical coefficient, in- 
troduced for a further convenience. The indices i,j label the faces of the "surface of initial 
data" discussed above. Since the commutation relations ([19|) are ultra-local (in the sense that 
the angle variables on different faces commute with each other), let us concentrate on the local 
Heisenberg algebra, 

H : [a,(]]=^h, (50) 

for a single lattice face (remind that the angles shown in Fig. [3] are related by (jl7p ). The 
map (|22p contains the quantities k,a,a*, defined in (jl9p . which now become operators. For 
definiteness, assume that the non-commuting factors in ()19p are ordered exactly as written. 
Then the definitions (|19p give 



a = q~h {U -U-^)-^ {UV -U-^V^^), (51) 

a* = {U -U-^)-^ {UV-'^ -U^^V), 

where the elements U and V generate the Weyl algebra, 

UV = qVU, C/ = e'", y = e*^, q = e'^^ . (52) 
The operators (|5ip obey the commutation relations of the gr-oscillator algebra. 



OSCg 



qa*a — q ^aa* = q — q ^, ka* = qa*k, ka = q ^ak 
P = g (1 - a* a) = q^^ (l-aa*) . 



(53) 



where the element k is assumed to be invertible. This algebra is, obviously, a quantum counter- 
part of the Poisson algebra (j24p . In the previous Section we have already mentioned the result 
of [21] that 

(i) the map (|22p is an automorphism of the tensor cube of the Poisson algebra (|24p (remind 
that the relation (1201) should be taken into account in (1221)1. 



In the same paper [21] it was also shown that 

(ii) there exists a quantum version of the map (j22p . which acts as an automorphism of the 
tensor cube of the g-oscillator algebra (j53p . The formulae (j22p for the quantum map stay 
exactly the same, but the relation (|20p should be replaced by either of the two relations 
on the second line of ([53l) . for instance, fc^ = q{l — a*a). In particular, ([23}) should be 
replaced with 

(A;^)2 = g(l-a*V) . (54) 
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For irreducible representations of the g-oscillator algebra (|53|) the formulae (|22p and 
uniquely determine the 7?.i23 as an internal automorphism, 

'R-i23{F) = Ri23 F Ru3, F e Osc, ® Oscg (g) Oscq . (55) 
It follows then from ()13p that the linear operator R satisfies the quantum tetrahedron equation 

-R123 -R145 ^246 ^356 = -R356 -^246 -^145 -^123 > (56) 



where the operators iii23, -R145, -R246 and -R356 act as (1551) in the three factors of the tensor 
product of the six algebras g-oscillator algebras (indicated by their subscripts) and act as the 
unit operator in the remaining three factors. 

4.2 Fock space representation model 

In this subsection we set ^ = 1 in (j49p . Assuming ?i > we have q = e^^ < 1. Define the Fock 
space representation of a single g-oscillator algebra ([53l) . 



F,: a|0)=0, |n) = -^=4= |0) , k\n) = q^+'/' \n) , n>0 (57) 
Then using (|22p . (j53p . (|54p and (|55p one can show that the the matrix elements of R are given 

by 



{ni,n2,n3\K\n^,n2,n3} - dn^+n^y^+n'^^^n^+n^y^+ni, Y [q^-q^]^^ {q^\q^)n2 {<P\<P)n 



X iziL^____ il- ,6Aq-^< O^Cl+n^) 2(l-n^+n3) 2 2(l+nl)^ 

(r;r)n' (r^^+"^^r)oo 



(58) 

(x; g2)„ = (1 _ a;)(l - <;2^) ... (1 - , (59) 

, . , 2 X (a;g^)n(fe;g^)n „ 

^^^^"'''^''^'^^ =S(<Z^;.^)„(c;,^)/ 



where 
and 



is the g-deformed Gauss hypergeometric series. Eq. (j55|) together with (p2|) and (jMj) lead to 
recurrence relations for the matrix elements of R. Such relations along with expressions for 
matrix elements of R for small values of the occupation numbers nj, were obtained in [21] (see 
Eqs. (30-32) therein). We have now solved those recurrence relations and obtained the explicit 
formula ([55|) . given abov^. This 3D iZ-matrix satisfies the constant tetrahedron equation ([56]) . 

n-^^,n^,n^ Yi'^^n^,n^ n'2,n'^,nQ n'^,n'^,n'Q / ^ "-S'^s'^e n2,n^,ng n-^^,n'^,n'^ n[,n'2,n'^ ^ ' 

n'.=0 n'.=0 

J J 



'Eq.JSHJ follows from Eqs. (29,30) of [21] where P^, (g^", g2^^ ^27) jg substituted with 

f 2(1-13'+-/) 2\ 

P ^2/3 27X _ (q -yQ )^ , , -213' 2(l-/3'+/3+7) 2(1-/3' +7) 2 2(1+q)x 

t^fs'Kl ,q ,<1 ) — I 2n+-Yi — 2^ — 2<PUg ,q ,q ,q ,q J 
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where the sum is taken over six indices n[,n'2,n'^,n'^,n'r^,ng and 

ld\%2^l = ("-!> "2, n3|i?|n'i, 712,^3) • (62) 

Note that Eq. (j6ip does not contain any spectral parameters. 

Now define a model of lattice field theory. On each edge of the cubic lattice place a fluctuating 
spin variable n, taking an infinite number of integer values n = 0, 1, 2, . . . 00. To each vertex of 
the lattice assign a local weight factor {ni^n2,n^\R |n']^, 71-2, "^3), given by ([58]) . depending on six 
spin variables placed on the edges surrounding the vertex, arranged as in FiglSJ The partition 



ni 




n2 



Figure 8: Graphical visualization of three-dimensional i?-matrix (j58[) and its matrix elements. The 

orientation of the axes is consistent with Fig. [2] 



function 

Zf = ^ ('^i) "^2, nsli? |ni, 712,71.3) , (63) 

(spins) (vertices) 

is defined as the sum over all spin configurations of the whole lattice where each configuration 
is counted with the weight equal to the product of the vertex weights over all lattice vertices. 
For definiteness we assume fixed boundary condition. The subscript "F" stands for the "Fock 
representation model" . 

In the quasi-classical limit 

q = ^1, h ^ 0, riih = — log ki = finite, (64) 

with large occupation numbers rii, n'^, such that ki and k'^ are kept finite, one obtains 

lim (7ii,7i2,7Z3|i? 171^,712,^3) = exp (- ^^^%^) ) (65) 

where the arguments of C^{k,k') are related by kik2 = k'^k^ and ^2^3 = A:2A;3. The function 
Cf{k,k') coincides with the classical Lagrangian (j45p provided the variables k,k' in (j64p are 
identified with those in (j34p . The sum (j63p in this limit can be replaced by an integral and 
estimated by the saddle point method, 

logZp = + 0(?i°), h^O, (66) 
n 

where S^'^''^ is the classical action (j46]l evaluated at its stationary point. 
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4.3 Modular double model 

In this subsection we set in P9]) 

^ = -i, h = Tih'^, (67) 

where 6 is a freeparameter, Re 6 ^ 0. Here it wih be more convenient to work with a shghtly 
modified versiorf^ of the map (|22|) . with the value e = — 1. It is worth noting that this modifica- 
tion does not afi'ect the bulk properties of the classical system (jUJ and leads only to boundary 
effects. In particular the equations of motion (j44|) remain unchanged. 

Consider a non-compact representation [44] of the (/-oscillator algebra (|53p in the space of 
functions /(fx) € L-^(IR) on the real line admitting an analytical continuation into an appropriate 
horizontal strip, containing the real axis in the complex a-plane (see [44] and [38] for further 
details) . Such representation essentially reduces to that of the Weyl algebra 

\Nq: kw = qwk, g = e*'^^', (68) 

realized as multiplication and shift operators 

k\a) = \e'^''^\a) , w\a) = \a-\b). (69) 

The generators a, a* in (f53l) are expressed as 

a= {l-qk^y/^w-\ a* = {l-q-^k^)^/^w. (70) 

As explained in [45] the representation (I69p is not, in general, irreducible. Therefore, the relation 
()55p alone does not unambiguously define the linear operator i?i23 in this case. Following the 
idea of [45] consider the modular dual of the algebra (f68l) , 

Wg-: kw = qwk, q = e-'^^'\ (71) 

acting in the same representation space 

k\a) = -\e'"'^'^ \a) , w \(t) = \a + \h-'^) . (72) 

We found that if the relation ()55p is complemented by its modular dual 

^123 {F) = iii23 F Rlis^ F e Oscg- ® Oscg (8) Osc,- , (73) 

then the pair of relations ([55]) and ([73]) determine the operator i?i23 uniqueljl^. The dual q- 
oscillator algebra OsCg is realized through the dual Weyl pair (|7ip and the relations 

a= (l-gP)l/2y)-l^ = {l-q-^P)^/^w. (74) 

The dual version of the map 7^i23 is defined by the same formulae ([22]) . where quantities kj, aj,aj, 
j = 1, 2, 3 are replaced by their "tilded" counterparts kj,dj, a*. The value of q does not, actually, 
enter the map ([22p . but needs to be taken into account in the relations between the generators 
of the g-oscillator algebra. Thus, the linear operator -R123 in this case simultaneously provides 
the two maps 7^123 and 7?.i23 (with e = —1). The explicit form of this operator is given below. 

^Note that the map (|22p with e — ±1 is a particular case a three-parameter map considered in [21]. 
®It is worth mentioning similar phenomena in the construction of the _R-matrix [46] for the modular double of 
the quantum group Uq{sl2) and the representation theory of Uq{sl2,V,) [47]. 
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Denote 

and define a special function 

'C1,C2 



6 + 5-1 



2*2 



/ Ci,C2 \ ^ /■ 
VC3,C4 J Jr 



(75) 



(76) 



wliere is the non-compact quantum dilogarithm [10] 



ip{z) = exp 



-2'\zx 



dx\ 

4 ^R+io sinh(j;6)sinh(x/6) x J 



(77) 



The values of ci, 02,03,04 are assumed to be such that poles of numerator in the integrand of 
(j76p lie above the real axis, while the zeroes of the denominator lie below the real axis. For 
other values of Cj the integral ()76p is defined by an analytic continuation. For Im6^ > the 
integral 2*2 can be evaluated by closing the integration contour in the upper half plane. The 
result reads 



2*2 



Cl,C2 
C3,C4 



Co 



jTT{co+ir]){ci—iri)—in{4:ri'^ +1) /12 



' C2 - Ci . 
7: 1- 



/ C3 - Cl \ / C4 - Cl 



X 2(^i( — qe''*^'^!-'^^)/''^ _ggi'(ci-C4)/b^ ~2 g7r(ci-C2)/f). ^2^ _ ~ g7r(c3+C4-ci -C2+2co)/f)^ 



(78) 



X2</'l(-'7e"('^-'^)^-ge^(^*-^l)^g2e-(^2-^l)^^7^-ge2-^^«) + (ci ^C2), 

where 2</'i is defined by (j60p . This formula is very convenient for numerical calculations. 
The kernel of i?-matrix serving the pair of the maps (j55p and (|73p is given by 



{(j\R\a') 



^'2^-2- 



, gi'f(o"i'^3+''?(°"l+'^3-'^2)) 



- c^a, -0"! + 0-3 
0-1 + 0-3, -o-i - 0-3 



0-2 



(79) 



It satisfies the constant tetrahedron equation ([56 



I da' ..da' R"'''"'^^"'' ^-p;.-^ 

i b '^l,<^2''^3 °'l''^4''^5 '^2''^4''^6 '^3''^5''^6 



i?'^^''^5.-6 ^-2.-4,-6 i^-l.-4--5 ^-1.- 
i D 0-3,(T5,(Tg 0-2,0-4, (Tg 0^,0-4,0-5 0-1,0, 



^11 ^" 

-.Oi ,0-2 ,03 



where 



n'^l>"2'<^i 

-''-0-1,02,03 



((Ti, 0-2, 0"3|i? |cr'i, 0-25 0-3) 

and the integrals are taken along the real line from —00 to +00 



(80) 



(81) 
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Similarly to (|63|) define a "modular double model" 



/■■■J n (^1^ 1^') n ' (82) 



(vertices) (edges) 



where the edge spins a now take continuous values on the real line. We assume fixed boundary 
conditions. Note that due to the presence of two delta- functions in ([75]) the edge spins are 
constrained by two relations 

CJi + O-g = cr'^ + (J3, + 0-3 = (T2 + (T3, (83) 

at each vertex of the lattice. 

Consider the quasi-classical limit of the i?-matrix (f79|) 



such that the variables 



6^0, \aj\, \a'A 00, (84) 



kj = ie^'^'^i , k'j = ie^'^'^J , i = 1, 2, 3, (85) 



are kept finite. Anticipating a connection with the classical formula ()45p let us identify the 
variables (j85p with those in (j34p . Further, for real b the squares k'j and k'^^ are negative (notice 
the factor i in front of the exponents in (|85p ). To accommodate this extra minus sign we will 
use a slightly modified form of (j35p . namely, 



. ^2 


^2\u, _ V2 


g2ia;o _ ""2 


g2i^3 _ ^2 


/ / ' 


1 ' 






_ V2 




„-2i^0 _ ^2 


g-2ia;3 _ ^2 




/ ' 




VlV\ 



(86) 



where the new variables Uj and a;^- obey the relations 

'^fe = - ('^o + ^^1 +^^2 + W3)/2, /c = 0,1,2,3. (87) 

Note also a useful identity 

uJi + ujj + u}'^+ uj[ = TT, (i, j, /c, Z) = perm(0, 1, 2, 3), (88) 

where (i, j, fc, /) is any permutation of (0, 1, 2, 3). If the cj-variables lie in the domain 

0<u;fc<7r, 0<a;fc<7r, /c = 0,l,2,3, (89) 

then in the quasi-classical limit (j84p the integral (j76p entering the formula (j79p can be estimated 
by the saddle point method 



(c7|i?|a') = exp|-^^^^i^ + 0(l)| , 6^0. 



(90) 



In practice, it is easier to calculate logarithmic derivatives of ()79p with respect variables kj and 
k'p rather than the expression (|79p itself. In this way one obtains 



- 1 ^ TT 

dCM{k,k') = -""^^ [ujd log sin ujj + ujjd log sin i^j) — — (^d log sin L02 + ■ (91) 



j=0 
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It follows then 

C\^{k,k') = £M(t^, 'w') + Legendre terms , (92) 

where the non-trivial part 

1 ^ 

Cuivy) = - Y,{H^3) + H^'j)) > (93) 

is expressed in terms of the Lobachevsky function [48], 

Jl{uj) = - log2sinxdx, < cj < vr , (94) 
^0 

while the Legendre terms read 

1 ^ TT 

Legendre terms = - [wj log 2 sin uj + uij log 2 sin Wj) — — (log 2 sin lj2 + log 2 sin 0^2) • (95) 

j=0 

Note, that in terms of the variables k,k',v,v' (connected with cj's by ([34]) and (I86jl ) the differ- 
ential of (p3]) reads 

d/:M(w,w') = (logy(ilogt;j -log-j^dlogw^- j , (96) 

while the Legendre terms take the form 

Legendre terms = (^og -4- log Vj — log ^ log v'j\ . (97) 

' j=i V I I / 

The formulae ([96j) and ([97]) are in a complete agreement with ([33]) and ([15]) . except the overall 
normalization and the replacement kj — > kj/\ and k'j kj/\ which is a gauge transformation 
resulting from a different of choice of the canonical variables. The total action 

4f(W)= E ^Mm),k'im)) (98) 

coincides (to within an overall factor and boundary terms) with an analytic continuation 
of (I46p into a regime where the variables fc| are negative and the variables Vj are unimodular. 
Substituting ([90]) into ([82]) one obtains 
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log Zm = + 0(1), b ^ (99) 

where S^P is the classical action ()98p evaluated at its stationary point. It is obvious from (|93p 
that this expression is real. 

Before concluding this section let us make a remark about the spatial symmetry properties 
of the modular model. It is convenient to introduce another kernel, 

{a\R\a') = {a\R \a') e^''K-'^2+'^3+'^i-4+'^3)/2 (iqo) 

differing from (j79p by a phase factor which introduces particular "external fields" but does not 
destroy the integrability of the model. The new kernel R obeys simple relations 

0-1,0-2,0-3 03,02, o^ -02,0^,-03 ■ 

which generate the whole cube symmetry group. 
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4.4 The "interaction-round-a-cube" formulation of the modular model 

An inconvenient feature of the modular representation model formulated above is that the 
edge spins are constrained by the relations (j83p . Here we re- formulate this model in terms of 
unconstrained corner spins, which also take continuous values on the real line. Fig. [9] shows an 
elementary cube of the lattice with the corner spins "a, b, c, d, e, f, g, /i" arranged in the same 
way as in [15]. The corresponding Boltzmann weight reads 




Figure 9: The arrangement of corner spins around a cube. 



3 

W{a\e,f,g\b,c,d\h;Ji,J2,T3;ai,32,33) = exp ai(cri + (Ti)r?(-1)^"^ 



1=1 



-^3) (102) 



where cTj, a[ are defined as 

(7i = g + f- a- b-Ji, a2 = e + g- a- c-J2, a3 = a + d- e- f + J3, 
a[ = c + d — e — h — Ti, a2 = b + d — f — h — T2 , a'^^ = g + h — b — c + T3 



(103) 



and satisfy the constrains (jSSp . The parameters Ti, T2, T3 are analogs of the spectral parameters 
in the Zamolodchikov model [14] and its generalization for an arbitrary number > 2 of spin 
states [16]. They are related by Tj = log[tan(6'j/2)], where 9j, j = 1, 2, 3, are the dihedral angles 
in the Zamolodchikov model (they are usually considered as angles of a spherical triangle). The 
parameters are similar to the linear angles in Zamolodchikov model (the side lengths of the 
above spherical triangle). However, unlike the Zamolodchikov model, the two sets of parameters 
Tj and in our case are totally independent. The parameters a, enter only the exponential 
prefactor in (I102p : they play the role of the external fields in the model. 

Consider spatial symmetry properties of the weight function (jl02p . The group of cube is 
generated by two elements: r (the reflection with respect to the diagonal plane (a, d, h, g) in 
Fig. [9|) and p (the 90° rotation around the vertical axis). They act on the spin variables and 
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parameters Tj, as follows 

r : (a|e/5r|6c(i|/i;Ti,T2,T3, 31,32, as) (al/e^lc^dl/i; T2, Ti, T3; 32, ai, 33), 

p : (a|e/5|6c(i|/i;Ti,T2, 73,31,32, 33) {g\cab\fhe\d; -Ti,T3, -T2;7r - 31,33,71 - 32). 

(104) 

On can show show that the weight function (|lU2p is invariant under these substitutions and, 
thus, possesses the full cube symmetry group. 




(105) 



Figure 10: Graphical representation for the tetrahedron equations for interaction-round-a-cube models. 
The weight function (jl02p satisfies the tetrahedron equation of the form [15] (see Fig. llOp. 

/ dz W{a4\c2, ci, C3|5i, 63, b2\z) VF'(ci I62, 03, ^i|c4, z, C6I&4) 

JR 

xW"{bi\z, C4, C3|a2, 63, h\c5) W"'{z\b2, b^, 63105, C2, cejai) = 

= / dzW'"{bi\ci,C4,,C3\a2,a4,a3\z)W"{ci\b2,as,a4\z,C2,CQ\ai) 
Jr 

X W'{ai\c2, z, C3|a2, 63, ai jcs) W{z\ai,a3, a2|c4, C5, C6I64), 
where the four sets of the spectral parameters are constrained as 

r, = j2, n' = -T3, vi = v„ rr = ri, t- = Ti, t- = -t; (loe) 

and 

3i — 3i - 3i, 3i —33 — 33, 32 — 32 — 32, 33 —33—32. 

The partition function is defined as 

Z'm= J ■■■ I llW{a\efg\bcd\h) J] da (108) 

cubes vertices 
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where the first product is taken over all elementary cubes and the integral is taken over all 
vertex spins of the lattice a,b,c,.... This model can be thought of as a non-compact analog 
of the generalized Zamolodchikov model [16], corresponding to an analytic continuation of the 
latter to complex values of the number of spin states N = —}p'. 

A precise connection between the two partition functions ([82]) and (|108p requires a detailed 
discussion of the boundary conditions in both models which is postponed to [38]. Here we only 
note that expression (I82p is particular case of (1108^ corresponding to a special choice of the field 
parameters 31,82,33 in (I102p . 

5 Conclusion 

In this paper we have exposed various connections between discrete differential geometry, sta- 
tistical mechanics and quantum field theory, displaying geometric origins of algebraic structures 
underlying integrability of quantum systems. 

We have shown that the 3D circular lattices are associated with an integrable discrete Hamil- 
tonian system and constructed two different quantizations of this system. The resulting 3D 
integrable models can be thought of as describing quantum fluctuations of the lattice geometry. 
The classical geometry of the 3D circular lattices arises as a stationary configuration giving the 
leading contribution to the partition function of the quantum system in the quasi-classical limit. 

We have also obtained two solutions of the tetrahedron equation, which naturally arise in 
our approach to the quantization of the circular lattices. One of these solutions is new. It has 
continuous spin variables taking values on the real line. This solution is connected with the 
modular double of the g-oscillator algebra and can be considered as a non-compact counterpart 
of the generalized Zamolodchikov model [16]. The other solution of the tetrahedron equation 
which arose here in the context of the circular lattices, was previously constructed in [21] by 
a purely algebraic approach. This solution possesses a remarkable property: it reproduces all 
two-dimensional solvable models related to finite-dimensional highest weight representations 
for all quantized affine algebras Uq{sln), n = 2,3, ... ,00 (the rank n coincides with the size 
of one dimension of the 3D lattice). Plausibly, a similar 3D interpretation, originating from 
other simple geometrical models, also exists for the trigonometric solutions of the Yang-Baxter 
equation, related with all other infinite series of quantized affine algebras [49, 50] and super- 
algebras [51]. Therefore, it might very well be that not only the phenomenon of quantum 
integrability but the quantized algebras themselves are deeply connected with geometry. 

Here we only stated our main results leaving detailed derivations to future publications 
[38]. There are many other questions we would also like to address there, in particular the 
geometric meaning of the Poisson algebra ()25p and connections of the 3D circular lattices with 
the 2D circle patterns [11] on the plane or the sphere. It would be interesting to understand 
underlying reasons of a "persistent" appearance of the g-oscillator algebra (j53p as a primary 
algebraic structure in many other important aspects of the theory of integrable systems, such 
as, for example, the construction of Baxter's Q-operators [52] and the calculation of correlation 
functions of the XXZ model [53]. It would also be interesting to explore connections of our 
results with the invariants of the 3D manifolds [54-56] , the link invariants [57-59] , quantization 
of the Techmueller space [60,61] and the representation theory of Uq{sl{2\M.) [62]. So, there are 
many interesting questions about the quantum integrability still remain unanswered, but one 
thing is getting is more and more clear: it is not just connected with geometry, it is geometry 
itself! (though the Quantum Geometry). 
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